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Vh ■ Grand canonical simulations are used to calculate adsorption isotherms of various 

P^ . classical gases on alkali metal and Mg surfaces. Ab initio adsorption potentials 
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, ^,. and Lennard-Jones gas-gas interactions are used. Depending on the system, the 
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resulting behavior can be nonwetting for all temperatures studied, complete wet- 
ting, or (in the intermediate case) exhibit a wetting transition. An unusual variety 



V^ . of wetting transitions at the triple point is found in the case of a specific adsorp- 

o, 

0^ , tion potential of intermediate strength. The general threshold for wetting near the 

"t^ . triple point is found to be close to that predicted with a heuristic model of Cheng 

et al. This same conclusion was drawn in a recent experimental and simulation 

Ch , study of Ar on CO2 by Mistura et al. These results imply that a dimensionless 

O 

wetting parameter w is useful for predicting whether wetting behavior is present 
^ , at and above the triple temperature. The nonwetting/ wetting crossover value 

. 5t , found here \s w c^ 3.3. 

PACS numbers: 64.70.Fx, 68.35.Rh, 68.45.Cd, 82.20.Wt 



I. INTRODUCTION 

Many adsorption systems exhibit either complete wetting below the triple point or triple 
point wetting transitions. These behaviors are recognized to arise from strongly attracting 
substrate potentials and/or small mismatch of substrate/adsorbate lattice constants [|l],3. 
Less attractive potentials exhibit incomplete or nonwetting behavior at the triple point. 
There is considerable interest in the nature of such weakly attractive situations near and 
above this temperature |^. General arguments and model calculations of Cahn |Q and of 
Ebner and Saam |^ in 1977 implied that systems which are nonwetting at the triple point 
ought to exhibit transitions to wetting at some temperature below the critical point. Experi- 
mental confirmation of this proposal was not found until appropriate systems were identified. 
In 1991, it was argued that alkali metals provide the weakest adsorption potentials of any 
surfaces for He atoms and therefore wetting transitions were expected for such substrates 
at tow temperature 0. This hypothesis was subsequently confirmed and the corresponding 
transition has since been studied in many laboratories 0^- More recently, similar phenom- 
ena have been predicted and/or seen for H2, Ne, and Hg on various surfaces ||T0|-[T5|. Key 
theoretical questions arising from these studies involve the reliability of both the assumed 
adsorption potentials and the statistical mechanical description of the transition, as well as 
the effects of surface heterogeneity [|TB| present in any real experiment. 

In the course of these investigations, Cheng et al. @ (CCST henceforth) posited a simple 
model which interprets the transitions in terms of a balance between the surface tension cost 
of producing a thicker film and the energy gain associated with the film's interaction with 
the surface, V(z) [0. There results an implicit equation for the wetting temperature (T^): 

{pi - p,)Iv = 27 (1) 



Iv = - dzV{z) (2) 

where pi and p^ are the number densities of the adsorbate liquid and vapor at coexistence, 



7 is the surface tension of the hquid, and Zmin is the equihbrium distance of the potential 
If one assumes a simple 3-9 form of the adsorption potential 

V = {AC^)/{27D^z'') - C/z^ (3) 

then the transition condition becomes 

(CD2)V3 = 3.337/(pz-p.) (4) 

CCST used this model to predict wetting transitions for a large number of specific systems. 
More recently, the same general approach was used by Chizmeshya et al. (CCZ) ||18|, who 
based their predictions on new ah initio adsorption potentials. The latter potentials differ 
from earlier potentials in several respects; especially important are differences in values 
of the coefficient C of the van der Waals asymptotic potential, resulting in significantly 
larger well depths D. There is indication that these revised potentials are more consistent 
with experimental data than are their predecessors. Particularly striking is the improved 
agreement for ^He contact angle p9|-pT| and wetting temperature on Cs, the most studied 



wetting transition to date J^-^]. 



These predictions have also been tested by recent experiments of Hess et al. |T^ and 
simulations of Ne on alkali metals and Mg [^. In this paper, we extend the simulation 
studies to a larger family of systems for which the hypothesis is relevant. We focus here on 
the following question: assuming that V(z) is known, does the CCST model correctly predict 
whether a specific adsorption system exhibits wetting at and above the triple temperature? 

The outline of this paper is the following. Section H summarizes our simulation method- 
ology and presents detailed results for several systems which yield fairly distinct adsorption 
behaviors. These include complete wetting for all T (Kr/Mg and Xe/Mg), an anomalous 
wetting transition near the triple point (Kr/PMg), a prewetting transition line (Ar/Mg) 
somewhat above the triple temperature, and prewetting transitions characterized by the 
formation of very thick films (for Ar, Kr and Xe on Li). 

Here PMg is our own nomenclature for "pseudo-magnesium", an ersatz surface which 



attracts Kr atoms about 7% more strongly than Mg itself [^. We also evaluate other 



hypothetical substrates in order to ascertain the triple point wetting threshold. In section |T|, 
we compare these findings with those obtained previously in simulations and/or experiments 



for inert gases on various surfaces [p^-p7| 



II. METHOD AND RESULTS 



The simulation technique has been described in detail previously [^, so only a brief 



description is given here. For this study, isotherms for various model systems were calculated 
using the Grand Canonical Monte Carlo (GCMC) technique |2^j2^. For each isotherm, 



6.8 million steps (each step being an attempted creation, deletion, or displacement of a 
molecule) were performed to reach equilibration and 4.5 million steps were performed in the 
data gathering phase. For some points near the transition regions, where the fiuctuations 
are larger, the number of equilibration steps was increased to 20 million and 12 million steps 
were performed for data gathering. 

As in our previous work, the model for the atom-atom interaction is a Lennard- Jones 
12-6 potential and the atom surface interaction is modeled using the potential described by 
Chizmeshya et al. ||l^. These potentials can be characterized by four parameters for each 



system studied. A well-depth egg and distance agg are needed to specify the adsorbate LJ 
potential and a well-depth D and the van der Waals constant C specify the atom-surface 
potential. The values of these parameters used in this study are listed in Table |I|. The 
parameters correspond to Ne, Kr, Ar, and Xe adsorbed on alkali and alkali earth metal 
surfaces. The surface labeled PMg (pseudo magnesium) was obtained when we attempted 
to simulate Xe on Li, but used the egg and agg values for Kr instead of Xe. This potential 
is similar to that of Kr on a Mg surface; however, the Kr/PMg distance of closest approach 
is 20% smaller and the well depth is about 7% weaker than estimated for Kr on Mg. While 
this potential does not correspond to a "real" system, we shall see that this small decrease 
in interaction strength significantly affects the wetting properties of the substrate/adsorbate 
system. 



In most cases the simulation cell had a height of h = 70 A, with the metal surface 
forming the boundary at one end and a hard wall surface at the other. In Xe/Li and Ar/Li 
simulations we increased the height to 140 A and 210 A. The surface unit cell was taken 
to be 10 (jgg X 10 agg givlug a nominal monolayer coverage of roughly 100 atoms. Periodic 
boundary conditions in these dimensions render the surface effectively infinite; we used a 
fairly large cutoff {5agg) to minimize long range corrections. 

For each system, we simulated adsorption isotherms over relevant temperature ranges 
from the triple point to the critical temperature. The results are shown in figures |1|-|^. Of 
the systems studied the least attractive surface of the surface interactions are for the Ar/Li, 
Kr/Li and Xe/Li cases. In all of these, we see evidence of prewetting transitions and, based 
on our results, we report wetting temperature. For Ar/Li we estimate T^ = 130 ± 4 K, for 
Kr/Li, T^ = 175 ± 4 K and T^ = 225 ± 4 K for Xe/Li. All T^ are reported in Table |. 

Despite the fact that these systems are examples of prewetting transitions, there are some 
qualitative differences from the case (say Ne/Mg) where the wetting transition appears close 
to the triple temperature. Figure |l] shows a series of isotherms for Xe/Li from T* = 1.0 to 
1.25. The isotherm at 221 K (the lowest T*) shows nonwetting behavior, manifested as slight 
adsorption below saturated vapor pressure. The higher temperature isotherms however, 
show prewetting transitions below saturated vapor pressure, from which we estimate Ty^. 
The nature of the prewetting jump can be further characterized by the results in figures 
0, 1^, and ^, all at T=254 K. Figure shows the dependence of the isotherm on the height 
h [= 70, 140, 210 A] of the simulation box. The rapid rise occurs at the same pressure, 
independent of box height. This is a clear indication that the transition is prewetting rather 
than capillary condensation. The detailed behavior is most easily understood by looking 
at the density profiles' dependences on h. Figure ^ shows these profiles at P = 22.9 atm. 
Note that the data for two largest boxes coincide but differ from those of the smallest box. 
Evidently, this regime requires quite large values of h in order to obtain reliable data; this is 
not surprising because the compressibility of the system is particularly high at this transition 



as the slope of the isotherms suggest. The same phenomenon is even more dramatic for the 
case of a 1% higher pressure. As seen in figure ^, one must go to even larger h ( > 200 
A) in order to derive accurate adsorption data even though the liquid regime occupies only 
the region z < 70 A. While this sensitivity to h and the details of the transition are quite 
interesting they are not the focus of this paper. Suffice it to say that the density jumps 
discontinuously near P = 23.0 atm, as seen in the different densities in figures § and §. 

We note one further point about these data which is that the film's density close to 
the surface is qualitatively similar to that on a relatively more attractive surface. There 
appears a quite large accumulation of adsorbate atoms in the region associated with the 
first layer. There is even a significant density maximum occurring where a dense second 
layer might occur. The latter is, however, preempted by the prewetting jump in the regime 
of temperature seen here. Above the prewetting critical point, in contrast, one does expect 
such smooth film growth. 

As the relative strength {D /egg) of the surface interaction increases, (Ar on Mg, Kr on 
PMg) the isotherms give evidence of wetting transitions. The results for Ar on Mg, (shown 
in figure ^) indicate prewetting transitions at about 10% above the triple point (T^^ = 81 
K). From these data, we estimate the wetting temperature to be 90 K for Ar/Mg with a 
prewetting critical temperature of 95 K. 

In figure ^, the results for Kr on PMg are shown. This system shows clear prewetting 
transitions below 128 K for all temperatures above the triple temperature (estimated to be 
116 K 10). Although this is a fictitious system, it represents one of the first simulations 
of triple point wetting. The only other simulation of wetting in this region is the case of 
Ne/Mg which we reported in ||22|| . 

Finally, as the surface interaction strength further increases, one begins to see continuous 
growth for all T between the triple temperature and the critical point. This is evident in 
figure 1^ which shows isotherms for Kr on Mg. The isotherms were simulated at temperatures 
near the estimated triple temperature and demonstrate complete wetting on the Mg surface 
by the adsorbate even at the lowest temperatures. The results for Xe on Mg (not shown) 
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are qualitatively the same. 

III. INTERPRETATION OF THE RESULTS 

To interpret our results we follow CCST, who evaluated their simple model in terms of 
the law of corresponding states. Then we note that equation ^ may be written in a universal 
form in terms of the reduced well depth D* and interaction range x, defined as follows: 

D* = D/e,, (5) 

X = {C/Df^la,, (6) 

The values of these parameters are reported in Table |. The CCST model may then be 
written in terms of a dimensionless wetting parameter 

w = xD* = 3.33-f/[{pi - Pv)(yggegg] (7) 

The triple point experimental data summarized in CCST yield a value of 3.6 for the right 
side of this equation. Hence the criterion is w = 3.6 for a wetting transition precisely at 
the triple point. The latter value, however, presupposes a 3-9 adsorption potential, which 
we have not used in the present simulations. If, instead, we compute the integral of the 
actual potential used in the calculations, we find a value which is 8 to 10% higher than the 
result obtained with the 3-9 potential. Taking this correction into account, the threshold 
for wetting at the triple point becomes somewhat lower for the realistic potential: 

w = 3.3 (8) 

This is shown in figure § as a dashed curve. Turning to the actual simulation results, we 
find that this prediction characterizes our data relatively well, insofar as the vicinity of the 
value 3.3 is a crossover regime of wetting behavior. For all gases exposed to the alkali metal 
surfaces the quantity w is less than 2.3, and the systems are nonwetting at the triple point. 
The case of Ar/Mg {w = 3.2) is one for which a wetting transition occurs about 10% above 



the triple temperature; similar behavior was found in our previous study of Ne/Mg {w= 
3.4). The case of Kr/PMg involves a similar value [w = 3.2), but we find anomalous triple 
point wetting transition behavior for that system. Finally, the cases of Kr/Mg [w = 3.4) 
and Xe/Mg {w = 3.8) both correspond to complete wetting. For comparison with a more 
typical adsorption system, the case of Ar/Au (w = 6.7) is one which exhibits wetting at and 
above the triple temperature i2^j27[ . 



We may add to this list of relevant systems the recently studied case of Ar/C02. Mistura 
et al. [^ found that it exhibits triple point wetting; this finding agreed with their density 
functional calculation using their theoretical potential. This potential would be roughly 
consistent with the CCST model in that the latter would predict a wetting transition at a 
temperature only 3% above the triple point temperature, 83.8 K |^ . 



While the criterion w = 3.3 successfully discriminates between the nonwetting and wet- 
ting data sets mentioned above, it is worthwhile to ascertain the validity of this criterion 
more accurately. Hence we have performed extensive simulations at the triple temperature, 
varying the well depth near the values x = 0.7, 0.9, and 1.4 in order to determine the wetting 
transition line. The results are shown in Figure §. One identifies this line as the point where 
the letter C is just above the letter W. Note that these transition values (dash-dot curve) lie 
somewhat above the w = 3.3 dashed line. This comparison indicates that the CCST model 
underestimates the well depth criterion for wetting. 

This finding is qualitatively consistent with results reported recently by Ancilotto and 



Toigo 1^. Using a density functional method, these workers found wetting transitions 
very close to the coexistence curve for several systems (e.g. Ne/Li and Na) for which 
our simulations find nonwetting behavior. The different findings of the simulations may 
be due to the uncertainty in identifying wetting behavior close to saturation, especially 
at T approaching T^. In these and other cases, the transitions occurred at (typically ~ 
10%) higher temperature than is predicted by the CCST model. These findings are also 
consistent with one's realization that the CCST model predicts wetting transitions for all 
adsorption circumstances, however weakly attractive, while a more reliable theory yields 
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drying transitions for ultraweakly adsorbing potentials p5[| . 

To summarize, we have found the CCST equation to be a semiquantitative guide to 
the occurrence of wetting transitions. It would be valuable to have additional, genuine 
experimental data to test this conclusion. We wish especially to encourage experimentalists 
to explore totally new systems, especially Ar, Kr, and Xe adsorption on alkali metals. There 
are several reasons for doing so. One is that such experiments will provide tests of the weak 
adsorption potentials for these systems, for which no experimental data exist. A second 
is that triple point wetting can involve interesting behavior at and just below the triple 
temperature, as has been found in many experiments [^^ . 



We have not explored either these systems, or analogous ones, below the triple tem- 
perature. In order to reliably describe the crystalline phase, as is necessary there, such 
simulations would require a more careful handling of the unit cell dimensions and shape 
than has been employed here. In view of the many relevant systems, such a study seems 
well justified for the future. 
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TABLES 
TABLE L The values of the LJ parameters, €gg and agg, for the gas-gas interactions and the 

values of the reduced well depth and width D* = D/cgg and x, as defined in Eq. (6), for various 

systems studied and correlated with wetting behavior, defined as follows: C = continuous growth, 

TP = anomalous wetting transition near or at the triple point (Ttr), W = wetting transition above 

the triple point, N = nonwetting for all T > Ttr, and D = drying. The Ar/C02 data were taken 



from Ref. [27|. The wetting temperatures estimated by the simulations are in parentheses. 



Gas / Surface 




Rb 


Na 


Li 


PMg 


Mg 


CO2 


Ne 


D* 


0.7 


1.1 


1.5 




2.8 




egg = 34 K 


X 


1.6 


1.4 


1.4 




1.2 




agg = 2.78 A 


behavior 


D 


N 


N 




TP(22K) 






xD* 


1.1 


1.5 


2.1 




3.4 




Ar 


D* 


1.1 


1.6 


2.1 




3.5 


3.8 


egg = 120 K 


X 


1.1 


1.1 


1.0 




0.9 


1.0 


agg = 3.41 A 


behavior 


N 


N 


W(130K) 




W(90K) 


C 
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1.2 


1.7 


2.1 




3.2 


3.8 


Kr 


D* 






2.3 


3.6 


3.8 




egg = 171 K 


X 






0.9 


0.9 


0.9 




agg = 3.60 A 


behavior 






W(175K) 


TP(116K) 


C 






xD* 






2.1 


3.2 


3.4 




Xe 


D* 






2.8 




5.45 




egg = 221 K 


X 






0.8 




0.7 




agg = 4.10 A 


behavior 
xD* 






W(225K) 
2.2 




C 
3.8 
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FIGURES 
FIG. 1. Solid curves are adsorption isotherms for Xe/Li at T=221, 232, 243, 254, 265, 276 K. 

The vertical dashed lines to the right of each curve indicate the saturated vapor pressures. 

FIG. 2. Adsorption isotherms for Xe/Li at T=254 K as a function of simulation cell height. 
The heights for the curves shown are 70 A(dashed), 140 A(dots), 210 A(full curve). The vertical 
dashed line to the right indicates the saturated vapor pressure. 

FIG. 3. Film density for Xe/Li at T=254 K and P=22.9 atm as a function of reduced distance. 
Curves correspond to varying heights of the simulation cell: h=70 A(dots), 140 A(dashed), 210 
A(dot-dashed), and 280 A(full curve). Note the latter three results coincide. 

FIG. 4. Film density for Xe/Li at T=254 K and P=23.1 atm as a function of reduced distance. 
Curves correspond to varying heights of the simulation cell: h=70 A(dots), 140 A(dashed), 210 
A(dot-dashed), and 280 A(full curve). Note the latter two results coincide. 

FIG. 5. Adsorption of Ar on Mg for T = 89, 90, 91, 92, 94, 95, 96, and 98 K (from left to 
right). The data indicate a wetting transition near 90 K and a prewetting critical point near 95 
K. The vertical dashed lines occur at saturated pressure for each temperature. 

FIG. 6. Adsorption of Kr on "pseudo-Mg", an artificial problem described in the text. The 
isotherms (from left) begin at 119 K and run to 128 K. There is a wetting transition close to the 
triple temperature, T^^ = 116 K. The prewetting critical temperature is 127 K. 

FIG. 7. Adsorption of Kr on Mg, showing continuous growth of the film for all T above the 
triple temperature. Curves correspond to T = 120, 125, 130, 135, 140, and 145 K, from left to 
right. 
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FIG. 8. Points represent systems shown in Table || or artificial systems studied as discussed 



in Section 01. Symbols are defined in the caption to this table. The threshold line based on the 
model discussed in the text, w = D*x = 3.3, is indicated by dashes. The dash-dot curve represents 
the wetting threshold derived from our simulations. 
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